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Abstract 

Recent experiments on URu2Si2 show that the low-pressure hidden order is non-magnetic but it 
breaks time reversal invariance. Restricting our attention to local order parameters of 5/ 2 shells, 
we find that the best candidate for hidden order is staggered order of either T@ or T xyz octupoles. 
Group theoretical arguments for the effect of symmetry-lowering perturbations (magnetic field, 
mechanical stress) predict behavior in good overall agreement with observations. We illustrate our 
general arguments on the example of a five-state crystal field model which differs in several details 
from models discussed in the literature. The general appearance of the mean field phase diagram 
agrees with the experimental results. In particular, we find that a) at zero magnetic field, there is 
a first-order phase boundary between octupolar order and large-moment antiferromagnetism with 
increasing hydrostatic pressure; b) arbitrarily weak uniaxial pressure induces staggered magnetic 
moments in the octupolar phase; and c) a new phase with different symmetry appears at large 
magnetic fields. 
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I. INTRODUCTION 



The nature of the so-called "hidden order" of the T < T ~ 17K phase of URu2Si2 
has long been debated Q|. Taking strictly on-site local order parameters only, U 4+ — > 5f 2 
shells can carry magnetic dipole, electric quadrupole, magnetic octupole, and even higher 
multipole order parameters. The full local symmetry is described by Q = T>ih®Qt where 
V±h is the tetragonal point group, and Q t = {E,T} the two-element group generated by 
the time reversal operator^] T. The classification of the twelve most obvious local order 
parameters is given in Table UJ Being expressed as Stevens equivalents, all order parameters 
are even under space inversion. The notation " g" and "ti" in Table refers to their parity 
under time reversal. 



TABLE I: Symmetry classification of the local order parameters |4j for B = (T>^h notations 



overline means symmetrization 
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About 30 years of work on one of the most intensively studied /-electron systems has not 
brought clarification: the order is still "hidden" As we are going to describe, theoretical 
progress has long been held up by the ambiguity of experimental findings on apparently 
heterogeneous samples. However, crucial recent experiments 0, Q] allow to infer what the 
equilibrium properties of ideal samples of URu2Si2 would be. 

From the earliest neutron scattering experiments 0, the issue has been complicated by 
the observation of apparent /-electron micromagnetism. Ascribing the magnetic moments 
to the bulk of the sample, the observations indicated two-sublattice Q = (0, 0, 1) antiferro- 
magnetism of U 5/-shell moments of O(0.01/xb) directed along the tetragonal fourfold axis 
z in the low-T (T < To) phase. Though the nominal value of the ordered moment m was 
two orders of magnitude lower than the paramagnetic moment, this seemed to conform to 
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the general idea that micromagnetism is the canonical behavior of /-electron systems on the 
borderline between the non-magnetic (heavy fermion) Kondo state and RKKY magnetism 
l|. According to this view, URu 2 Si 2 might have been put in the same class as UPt 3 or 
CeAl 3 Q- 

Many previous ideas about URu2Si2 were based on the assumption that antiferro mag- 
netism with micro-moments is a static phenomenon, and it is an intrinsic feature of the 
T < To phase. Since the ordering of small moments could not account for the large thermal 
anomalies at the 17K transition, it was assumed that the staggered dipole moment m is 
a secondary order parameter, being induced by the primary ordering of an un-identified 
full-amplitude order parameter if) (the hidden order). This would require that antiferro- 
magnetism has the same symmetry as the hidden order, i.e., if) should break time reversal 
invariance and share the spatial character of m under the symmetry classification according 
to the tetragonal point group T> ih . With these assumptions, the Landau free energy 
functional would contain a term —mif), generating m 7^ whenever the primary ip 7^ 0. 
This is a scenario which we are going to discard, for the reasons given below, and further in 
Sec. ED 

The intimate connection between hidden order and micromagnetism looked always some- 
what suspicious because the variability (0.017 — 0.04)/xb of the antiferromagnetic moment 
was too large to be associated with nominally good-quality samples, and because the onset 
of micromagnetism did not exactly coincide with To. Susceptibility and NMR under pressure 
give an insight: though the sample-averaged sublattice magnetization grows with pressure, 
it seems to arise from the increase in the number of magnetic sites, not from changing the 
magnetic moment at a given site . This points to the possibility that the apparent micro- 
magnetism is an attribute of heterogeneous samples, and should be understood as ordinary 
antiferromagnetism of a small (~ 1%) volume fraction in samples which for some reason 

n 

always include a minority phase [ L 13J . 

The argument was clinched by high-pressure /xSR experiments: hidden order is non- 
magnetic, and antiferromagnetism of at least O(0.1/ib) ionic moments appears at a first-order 
transition at p tr « 0.6GPa (Ref. la). There are two thermodynamic phases, a non-magnetic 
phase with (if)) 7^ and (m) = 0, and the antiferromagnetic phase with (m) 7^ and (if>) = 0. 
At ambient pressure, the magnetic ((m) 7^ 0) phase is slightly less stable than the phase with 
hidden order ^jj. However, large-amplitude antiferromagnetism is stabilized at hydrostatic 
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pressures p > 0.6GPa following a first order non-magnetic-to-magnetic transition. In a range 
of low hydrostatic pressures, the nature of the low-temperature phase remains the same as 
in ambient conditions: if) ^ and m = 0. The situation is, of course, different if we apply 
fields which lower the symmetry of the system: magnetic field B, or uniaxial stress a. 

In the following Sections, we discuss the effect of uniaxial stress, and of magnetic field, 
on the ordered phases of URu2Si2- We will deduce that the low-pressure zero-field order 
must be staggered octupolar order of either B\ u or B 2u octupoles (Sec. |Hj). The overall 
appearance of the temperature - magnetic field phase diagram will be explained (Sec. IIII|) . 
Finally, the general arguments will be illustrated by the results obtained from a new crystal 
field model (Sec.HVJ). 



II. OCTUPOLAR ORDER 

In this Section, we argue that the experimental evidence presented in Refs. |2l unam- 
biguously shows that the "hidden order" of URu2Si2 is alternating octupolar order with 
3 = (0,0, 1). Here we restrict our attention to the strictly local (on-site) order parameters 
15| listed in Table |lj Two-site quadrupole-spin and three-site spin-spin correlators could 
appear in the same symmetry class as on-site octupoles Q| ; the present argument does not 
differentiate these cases. 

Let us recall the hydrostatic pressure - temperature phase diagram obtained from high- 
pressure /iSR experiments p (a phase diagram of the same shape results from our mean 
field theory, see Fig. El left). Hidden order (ip) ^ is the attribute of the low-pressure, low- 
temperature phase (p < p tT , T < To(p)). Though all samples show some micro magnetism, it 
can be safely concluded that this is an extrinsic effect and in a perfect sample, hidden order 
should be non-magnetic. Antiferromagnetism of at least O(0.1/xb) ionic moments appears 
at a first-order transition at p tT ~ 0.6GPa . At p < p tr the hidden order onset temperature 
T (p), and at p > p tr the Neel temperature T N (p) are critical temperatures for the hidden- 
order-to-paramagnetic, and the antiferromagnetic-to-paramagnetic transitions, respectively; 
these meet the first order phase boundary at a bicritical point. It follows that m and if) are of 
different symmetry, and the Landau free energy cannot contain a term mip. The symmetry 
of if) must be in any case different from A 2u (Q = (0, 0, 1)). 

Regarding the absence of magnetism as an established fact, we can also exclude E u 
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(magnetic moments perpendicular to the tetragonal z-axis). The remaining choices for the 
order parameter ip are: quadrupolar B 2g , or E g ), octupolar (B\ u or B 2u ), hexadecapole 
(A 2g ), or triakontadipole (A lu ) (see Table HJ). Quadrupole and hexadecapole moments are 
time reversal invariant, while octupoles and triakontadipoles change sign under time reversal. 

An important recent experiment allows to decide the time reversal character of the hidden 
order. Yokoyama et al |7[ carried out magnetic neutron scattering measurement in the 
presence of uniaxial stress applied to a single crystal sample either along, or perpendicular 
to the tetragonal axis. Stress a || (001) does not produce significant change in magnetic 
moments. However, for stress o _L (001) the staggered moment increases approximately 
linearly, reaching ~ 0.25/iB at a = 0.4GPa. In contrast to hydrostatic pressure, no threshold 
value is needed to induce a magnetic moment; it appears as soon as the stress a is finite. 

Mechanical stress is time reversal invariant, thus it can produce magnetic moments only 
from an underlying state which itself breaks time reversal invariance. This limits the choice 
of hidden order to B\ u or B 2u (octupolar), or A\ u (triakontadipoles). We emphasize that the 
choice of octupolar order is essentially different from the previously assumed quadrupolar 
order 0,^,0] which does not break time reversal invariance. Additional evidence in favor 



of the time reversal invariance breaking character of ip comes from NMR measurements) 

We will show that the properties of URu2Si2 can be described well with the assumption of 
octupolar order. This would make URu2Si2 the third well-argued case of primary octu pol ar 
order in an /-electron system (the first two cases being NpC>2 2l| and Cei-xLa^Bg 
Within the limits of our argument, either B\ u or B 2u would reproduce the basic effect of 
stress-induced large-amplitude antiferromagnetism. On the other hand, we rule out A\ u 
triakontadipoles as order parameters. 

First, we consider stress applied in the (100) direction, a || (100) lowers the symmetry to 
orthorombic V 2h (see Appendix). Under V 2h , Af^ — > B°^ h and B^ 1 — > -B°„ th , so the order 
parameters Tf and J z become mixed (Table ITVj) . A state with spontaneous Tjr octupolar 
order carries J z magnetic dipole moments as well j^, accounting for the observations [z]. 

An alternative way to derive this is by inspecting the relevant terms of the Landau 
potential for the undistorted tetragonal phase (the operators in the equations below have 
the meaning given in Table |1J). Choosing B 2u octupolar order parameter, consider the mixed 
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FIG. 1: Stress- induced magnetic moment in the octupolar phase, based on the crystal field model 
described in Sec. IW1 Thick line: (M z ) staggered magnetization, thin line: (Tf) octupolar moment, 
as a function of the uniaxial pressure a || (100) (a in arbitrary units). 

third order invariant 

X{A 2u ®B lg ®B 2u ) = Cl J 2 (0)T/(Q)O 2 2 (-Q) 

+c 2 J 2 (Q)T/(-Q)O 2 2 (0). (1) 

Generally, c\ ^ and c 2 ^ 0. For our present purpose, the second term matters. A 
uniform stress a || (100) induces uniform (q = 0) 0\ quadrupole density pj] which cou- 
ples the staggered (q = Q) B 2u octupole order parameter to the J z dipole density with 
q = — Q, i.e., the same spatial modulation. Neutron scattering shows that stress-induced 
antiferromagnetism has the same simple two-sublattice structure with Q = (0, 0, 1) that was 
previously ascribed to micromagnetism, thus the periodicity of the hidden octupolar order 
must be the same. 

The actual stress-dependence of antiferromagnetic polarization depends on microscopic 
details. Fig. [T] illustrates the general behavior for a crystal field model which we describe in 
detail later. 

Stress applied along the z-axis induces 0\ which transforms according to the identity 
representation Ai g , thus it does not appear in the invariants, and it is not predicted to 
induce magnetism. This is in qualitative accordance with the observation that for o || (001) 
the induced moments are an order of magnitude smaller than for a || (100). We believe that 
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the fact that these moments are not exactly zero, is due to non-ideality of the sample, as 
micromagnetism itself is. 

The situation is less clear with varying the direction of stress in the a _L (001) plane. 
Experiments find that the stress-induced antiferromagnetic moment is essentially the same 
for <t||(110) as for er||(100) |7j. Taken in itself, stress-induced antiferromagnetism would be 
as easy to understand for <t||(110) as it was for <r||(100). Namely, the invariant expansion of 
the Landau potential contains also 

l(A 2u ®B 2g ®B lu ) = c 3 J z (0)T xyz (Q)O xy (-Q) 

+c 4 J z (Q)T xyz (-Q)O xy (0) . (2) 

cr || (110) induces uniform O xy quadrupolar polarization. Assuming that the hidden (octupo- 
lar) order is T xyz (— Q), it is coupled to J Z (Q), the same kind of antiferromagnetism as we 
found before. An alternative way to arrive at the same result is by observing that <t||(110) 
lowers the symmetry to orthorombic, under which T xyz and J z belong to the same irrep 
(Appendix, Table HVj) . 

We have to emphasize, though, that assuming a homogeneous system, either we have 
an explanation for the effect at cr || ( 100) (with Tf octupolar order), or for er||(110) (with 
T xyz octupolar order), but not for both. Under tetragonal symmetry, T~P and T xyz belong 
to different irreps, and therefore these orders cannot coexist. At the level of our present 
argument, the problem cannot be resolved. We believe that it is not merely a difficulty with 
our model but it points to a genuine feature of URu 2 Si 2 . We speculate that the T@ and T xyz 
orders are sufficiently near in energy, and so samples tend to contain domains of both. 

We note that the A\ u triakontadipole J x J y J z (J x — J y ) (see Table HJ) would not give rise 
to stress-induced magnetism and is therefore not a suitable choice as order parameter. 

It is worth pointing it out that our present scenario offers an explanation why micro- 
magnetism is always present. This may seem paradoxical since if it were connected with a 
minority phase only, it would be reasonable to expect that some preparation techniques give 
single-phase samples, i.e., completely non-magnetic ones. However, we ascribe antiferromag- 
netism also to the polarization of the primary octupolar phase in a stress field. It can be 
assumed that the environment of impurities and crystal defects always contains regions with 
the local stress oriented perpendicularly to the tetragonal main axis, thus there is always 
some local antiferromagnetism. 
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III. MAGNETIC FIELD 



There have been extensive studies of the effect of an external magnetic field on the phase 
digram of URu 2 Si 2 QQ. The Sy ste m is lively h_e to holds applied in the X -y 
plane, while fields B || z have substantial effect: hidden order can be suppressed completely 
with i?cr,i = 34. 7T. The phase boundary in the B-T plane is a critical line, thus hidden 
order (or its suitable modification) breaks a symmetry also at B ^ 0. At somewhat higher 
fields, an ordered phase appears in the field range B CTt 2 = 35. 8T < B < £? cr 3 = 38. 8T. One 
possibility is that it is the re-entrance of the B < B cr i hidden order; however, we are going 
to argue that the high-field order has different symmetry than the low-field order. 

The difference between the previously suggested quadrupolar order [itJ , and our present 
suggestion of octupolar order, is sharp at B = (Ref. |2?). However, a5^0 magnetic field 
mixes order parameters which are of different parity under time reversal (Table ITT|) . The 
reason is that switching on a field B || z lowers the point group symmetry from T>^f l ®Q t to 
an 8-element group isomorphic (but not identical) to C$ v (Ref. I4). 

Switching on a field B || z, geometrical symmetry is lowered from T>^ to C^. However, 
the relevant symmetry is not purely geometrical. Though taken in itself, reflection in the xz 
plane a v , x is not a symmetry operation (it changes the sign of the field), combining it with 
time reversal T gives the symmetry operation Ta VjX . The same holds for all vertical mirror 
planes, and C 2 1 2 axes, thus the full symmetry group consists of eight unitary and eight 
non-unitary symmetry operations [3] 

Q{B Z ) = C 4h + f(7 VjX C 4h . (3) 

We may resort to a simpler description observing that 

Q = C 4 + fa VyX C A (4) 

is an important subgroup of Q(B Z ), and we can base a symmetry classification on it. The 
multiplication table of Q is the same as that of C^ v , and therefore the irreps can be given 
similar labels. It is in this indirect sense that the symmetry in the presence of a field B || z 
can be regarded as C± v (a convention used in Ref. |5j). The symmetry classification of the 
local order parameters valid in B || z is given in Table ITT1 The results make it explicit that 
the magnetic field mixes dipoles with quadrupoles, quadrupoles with certain octupoles, etc. 



8 



TABLE II: Symmetry classification of the lowest rank local order parameters for B || z (notations 
as for C± v |2() 



Symmetry 


basis operators 


Ai 


1, Jz 


A 2 


Jx<Jy(Jx t ^y)i Jx<Jy<Jz{ Jx Jy ) 


Bi 


Ol Tf 


B 2 


O T 

^xyi ■'■xyz 


E 


{JxiJy\i {Qxzi@yz\ 



In a field B || z, there can exist ordered phases with four different local symmetries: A 2 , 
Bi, B 2 , and E. The zero-field I^w-type Tf octupolar order evolves into the .Bi-type Tf- 
0\ mixed octupolar-quadrupolar order. If the octupolar order is staggered, it mixes with 
similarly staggered quadrupolar order: this follows from the first line of Eqn. ((TJ The 
character of the low-field phase is indicated in the ground state phase diagram in Figure El 
(all numerical results are derived from a crystal field model described in Sec. IIVI but the 
validity of our general arguments is not restricted to that particular model). The gradual 
suppression of octupolar order under field applied in a high-symmetry direction is a well- 
known phenomenon; a similar result was derived for f 3 ions in Ref. |29j. In our calculation, 
the octupolar phase is suppressed at B crl « 34. 7T (Figure EJ). 

Although we are not familiar with experimental results for the combined effect of hydro- 
static pressure and magnetic field, it should follow from our scheme that a critical surface 
is bounding the phase with staggered B\ octupolar-quadrupolar order until at sufficiently 
high pressures, the critical surface terminates by a bicritical line. The high-pressure low- 
field phase has alternating J z order like in the zero-field case. Hydrostatic pressure does not 
change the symmetry of the system, but it can change the numerical values of the coeffi- 
cients in the expansion of the Landau free energy in terms of invariants. Therefore, generally 
speaking, we expect continuity with the results found for p = latm up to a threshold value 
of the pressure where a first order transition to a phase with different symmetry may take 
place. 

Let us return to the case of B || z field effects at ambient pressure. The story of the 
gradual suppression of the B\ octupolar-quadrupolar phase is closed by itself; it might have 
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FIG. 2: The high-field part of the T = phase diagram of the multipolar model (B in units of T 
(Tesla)). Vertical axis: (Tf) for the low-field phase, and (O zx ) for the high-field phase. The field- 
induced mixing of the order parameters is shown within the shaded areas. The overall appearance 
of the T-B (inset, T in units of K) phase diagram is very similar. (The critical temperature of the 
E phase is scaled up 3-fold). 

happened that there is only one ordered phase, surrounded on all sides by the disordered (Ai) 
phase. However, as shown in Table ITT) there are order parameters of different (A2, B2, E) 
symmetries; it depends on microscopic details whether such orders are induced by sufficiently 
high fields. If they are, they cannot coexist with Bi, so the corresponding domains in the 
B z -T plane must be either disjoint from the B±, or, if they are pressed against each other, 
separated by a first order phase boundary. Our Figure |2] illustrates the former case, where 
an E phase with mixed quadrupolar-dipolar order (see Table |H]) is separated from the low- 
field B\ phase by a narrow stretch of the disordered phase. We observe that this model 
result close resemblance to the phase diagram determined by high-field experiments 

I25L [26 1 . We note that high-field transport experiments add more phase boundaries to those 



determined by static experiments However, it is often found that transport anomalies 
delineate regions which, while showing interesting differences in the dominant conduction 
mechanism, still belong to the same thermodynamic phase. Therefore, we take the view 
that the boundaries shown in Figure El are the most robust features of the phase diagram, 
and the first step should be identifying the nature of these. 
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FIG. 3: In the low-field regime, the minimum model consists of three singlets. Left: The field- 
dependence of the levels. Right: relevant multipole matrix elements. 

An interesting possibility to recover a phase diagram of the same shape would be to 
identify the high-field phase as the "re-entrance" of the low-field Bi phase. This possibility 
was suggested in Ref. [jj. However, our present model study does not predict re-entrance. 



IV. CRYSTAL FIELD MODEL 



The previous arguments were based on a symmetry classification of the order parameters, 
and the conclusions are independent of the details of the microscopic models that allow the 
emergence of the ordered phases (in particular, zero-field octupolar order) which we pos- 
tulate. However, many physical properties (foremost the temperature dependence of the 
susceptibility, but also the specific heat) were fitted with apparent success by making differ- 
ent assumptions about the nature of the zero-field hidden order (either quadrupolar order 
of 5/ 2 shells [n, IS, 3(3], or non-conventional density waves |m|)- Therefore it is important 
to show that our work is not in conflict with findings for which alternative explanations had 
been suggested but offers fits to the results of standard measurements, which are at least 
comparable, and in some cases better, than previous results. 

Here we assume that equilibrium phases other than the superconducting phase, can be de- 
scribed in terms of localized /-electrons, with stable 5/ 2 shells. We note that for many other 
interesting /-electron systems (e.g., CeB 6 and Pr-filled skutterudites) the localized-electron 
description of multipolar ordering works well, in spite of the fact that for certain physical 
quantities, consideration of the itinerant aspects of /-electron behavior is indispensable. 

It is generally agreed j^] that the crystal field ground state is a singlet, and that the 
salient feature of the level scheme is three low-lying singlets. Three singlets are sufficient to 
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TABLE III: Tetragonal crystal field states used in the model. 



state 


form 


symmetry 


energy [K] 


1*2) 


1/V2(|4) - |-4» 


A 2 


100 


\d±) 


a\±3) - Vl - a 2 |=Fl> 


E 


51 


\u) 


1/V2(\2) - |-2)) 


B 2 


45 


N 


6(|4) + |-4)) + Vl-26 2 |0) 


A x 









112^ 
\ A . ~> 








1 t 4 > 


o yz ,o 2S 

J x ' ^ y 


1 



FIG. 4: Left: The magnetic field dependence of the single-ion levels in the extended five-state 
model used up to high values of T and B. Right: Additional multipole matrix elements due to the 
addition of the doublet state to the crystal field levels, which are relevant for the high-field phase. 



account for low-energy phenomena. It is found that further two states have to be taken into 
account to get a satisfactory fit for the susceptibility up to room temperature. We note that 
the nature of the high-field ordered phase has not been discussed in previous crystal field 
theories. 

The backbone of our crystal field model is the inclusion of the same three singlets as in the 
works of Santini and coworkers 12, QQ, but in different order (Table HH Fig. El (left)). 
The ground state is the \ti) singlet, and fa) an excited state lying at A 2 = 100K. 1^) and 
^2) are connected by a matrix element of J z , as observed by neutron inelastic scattering 0. 
The lower-lying singlet ^4) is connected to the ground state by an octupolar matrix element: 
this feature allows the existence of induced octupolar order as the strongest instability of the 
system (Fig. [21 right). We remark that while other level schemes may also allow octupolar 
order if one assumes a stronger octupole-octupole interaction, our assumption seems most 
economical. 

Finally, as in previous schemes, at least two further states are needed to fit magnetization 
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data up to 300K. We found it useful to insert one of the doublets (\d±)). This is an alternative 

nn 

to models with five singlets jl 71 . 130| . As we are going to see, fits to standard macroscopic 
measurements are no worse in our scheme than in previous ones. However, our scheme has 
the advantage that it accounts for the high-field observations. We show the field-dependence 
(B || z) of the crystal field levels in Fig. 0] The salient feature is the crossing of the (mildly 
field- dependent) singlet ground state with one of the levels derived from the splitting of the 
doublet at a field strength lying between the critical fields B cr ^ and _B cr 3. The crossing levels 
are connected by matrix elements of E operators ( Table ITT]). Consequently, we find a high- 
field E phase where { J x , J y }-type transverse dipolar order is mixed with {O zx ,O yz }-tjpe 
quadrupolar order (see Figure |2J). 

Commenting on differences between our crystal field scheme (Table 11111 where we use 
a = 0.98, b = 0.22) and previously suggested ones, we note that unambiguous determination 
is very difficult even if an intense experimental effort is undertaken, as in the recent case of 
Pr-filled skutterudites. By and large we agree with Nagano and Igarashi [3] , who argue that 
the crystal field potential of URu2Si2 is not known in sufficient detail yet. We complied with 
constraints which appear well-founded, as e.g. the neutron scattering evidence by Broholm 
et al 8], but otherwise we adjusted the model to get low-field octupolar order for which 
we found model-independent arguments. Level positions were adjusted to get good overall 
agreement with observations but we did not attempt to fine-tune the model, neither did we 
check for alternative schemes with less straightforward parametrization. 

We use the mean field decoupled hamiltonian 

Hmf = A 1 |t 4 )(t 4 | + A 2 |t 2 )(t 2 | + A 3 K)(d a \ 

a=+- 

-gfi B BJ z + A oct (Tf) Tf - A quad (O zx ) O zx (5) 



where g = 4/5, and the octupolar mean field coupling constant A oct is meant to include the 
effective coordination number; similarly for the quadrupolar coupling constant A qua( j. We 
assume alternating octupolar order, and uniform O zx order; the result would be the same 
if the high-field quadrupolar order is also alternating. We do not introduce independent 0\ 
or { J x , J y } couplings, nevertheless (0%) 7^ in the Bi phase, and (J x ) 7^ in the E phase. 

At B — 0, the only non-vanishing octupolar matrix element is C = (ti|7/|t4) ~ 8.8. 
Octupolar order is driven by the large C: assuming A oct = 0.336K we get the critical 
temperature T (B = 0) = 17.2K for 7/-type antiferro-octupolar order. Using a similar 
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estimate, we find A„ c p ~ 0.2K for NpC>2 which orders at 25K (Ref. l2ll ). The order-of- 
magnitude correspondence between two documented cases of octupolar order shows that 
our present estimate of the octupolar coupling strength is not unreasonable. 

(JHJ) was solved for all temperatures and fields B z . We find that the octupolar phase is 
bounded by a critical line of familiar shape (Fig. [21 inset), which has its maximum T = 17. 2K 
at B z = 0, and drops to zero at -B C r,i = 34. 7T. The transition remains second order through- 
out; we did not hit upon a tricritical point, though we are aware of no reason of why it 
should not have appeared. 

Similarly, the ground-state amplitude of the octupolar order is a monotonically decreasing 
function of B z (Fig.|2J). 

The restricted model with three singlets (Fig. ^ offers two basic choices. In the absence 
of symmetry-lowering fields, the Landau expansion of the free energy in terms of the order 
parameters is 

T = ao(T,p)(T^ 2 + Po(T, P )(T z ^ 4 + a M (T,p)(J z ) 2 

+P M (T,p)(J z ) 4 + ... (6) 

Note that because of the tetragonal symmetry, the free energy expansion does not contain 
the term (T^)(J Z ). It follows that the possible ordered phases can be (A) (Tf) ^ and 
(J z ) = or (B) (T£) = and (J z ) ^ 0. This is in agreement with the experimental finding 
p that (A) is the low-pressure phase, and (B) is the high-pressure phase, and they are 
separated by a first-order boundary. 

This canonical case (qualitatively agreeing with the schematic phase diagram shown in 
Ref. s| is illustrated in Fig. 03 (left). It was derived from Eqn. (JHJ) using an ad hoc model 
assum ptio n about the pressure dependence of the crystal field splittings Ai and A 2 (Fig. 
right) [34j. The shape of the phase boundaries, and in particular the slope of the first order 
line, could be fine-tuned by adjusting the pressure dependence of the crystal field parameters, 
but the overall appearance of the phase diagram: two critical lines meeting at a bicritical 
point, which is also the end-point of a first-order boundary, is generic. 

The stress dependence of the induced antiferromagnetic moment (Fig. was determined 
in a similar calculation, adding the term — aO\ to the Landau potential, and solving the self- 
consistency equations for (J z ) and (Tf). a in this calculation has the character of uniaxial 
stress, but an additional set of experimental data would be needed to determine its absolute 
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FIG. 5: Left: The pressure dependence of the critical temperature of the octupolar and the dipolar 
antiferromagnetic phases (Tq(p) and Tjv(p), respectively). The first-order boundary between the 
two ordered phases is an interpolation through the calculated points. Right: model assumption 
about the pressure dependence of the crystal field splittings Ai and A2. 
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FIG. 6: Linear susceptibility per site (in //b/T) on extended temperature scale (left), and in the 
vicinity of the octupolar transition (right). The dashed line gives the single- ion result. 

scale. 

Next, we consider the results of some standard low-field measurements. This was not 
the primary purpose of our work but rather serves as a check. The quadrupolar model ^| 
obtained a reasonably good fit for the temperature dependence of the linear and non-linear 
susceptibility in a range of temperature, and we have to prove that our model yields a 
comparably good description on a completely different microscopic basis. 

The octupolar transition shows up as a discontinuity of the linear susceptibility (Fig. EJ). 
This character of the xi anomaly is expected from general arguments (3^. While the low- 
temperature behavior, including the regime around T , is satisfactorily described by the 
three-state model (Fig. 0J), fitting the susceptibility up to room temperature (Fig. left) 
requires the five-state model (Fig. HJ). O ne of the hallmarks of the hidden-order transition 
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FIG. 7: The temperature dependence of the nonlinear susceptibility X3 m the vicinity of the 
octupolar transition. The dashed line is an interpolation through the calculated points 
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FIG. 8: The magnetization curve at T = (M in units of /ib)- 

of URu2Si2 is the strong jump of the non-linear susceptibility X3 (Ref. 0). The shape of 
the calculated anomaly (Fig. EJ) corresponds rather well to the experimental result. 

Next we discuss the high- field behavior at B > B cr \, and interpret the disjoint high- 



field phase observed in experiments [25. 26] mixed quadrupolar-dipolar phase (the E 
phase in Fig. |2J). We exploit the field dependence of the ionic levels in the five- level model 
( Table ITTT1 Fig. HJ). The single-ion levels t\ and gL would cross at -E> cr0 ss = 37.3T. Since 
\t\) and \d-) are connected by E operators including O zx (see Fig. HI right), a range of 
fields centered on -B C ross is certain to favour {O zx , O yz \ quadrupolar order, and simultaneous 
{Jx, J y } dipolar order. We chose a weak quadrupolar interaction A quac j = 0.054K in Eqn. 
this gives quadrupolar order between the critical fields B cr 2 = 35. 8T and B ct 3 = 38. 8T. The 
amplitude of quadrupolar order is not small (Fig. |2J) but the ordering temperature is low 
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(~ IK) because the coupling is weak. The E phase shows up as the steep part of the 
magnetization curve in Fig. |HJ For A qua( j = we would have a jump-like metamagnetic 
transition at B = B cross . 

We are aware of an unsatisfactory feature of the calculated magnetization curve. Though 
it is clear that our theory involves three critical fields: B CT i, B cv2l and B cv 3 , at the lowest 
of these the anomaly is so weak that it does not show up on the scale of Fig. |HJ We get a 
single-step metamagnetic transition distributed over the width of the high-field quadrupolar 
phase. The overall height of the step is right, 
of the transition observed by Sugiyama et al 

V. DISCUSSION AND CONCLUSION 



3ut we do not recover the three-step structure 

3- 



There have been many attempts to explain the non-superconducting phases of URu2Si2. 
Though the behavior of /-electrons in this system certainly has itinerant aspects, or perhaps 
URu2Si2 is on the verge of a localized-to- itinerant transition, arguing on the basis of a simple 
localized electron model can lead to useful results. Namely, crystal field theory conforms 
to a general symmetry classification of the equilibrium phases, which is expected to apply 
to a wider range of models, including suitably defined Kondo lattice, or Anderson lattice, 
models. Our main interest lies in cross-effects like the mixing of order parameters in the 
presence of external magnetic field, or mechanical stress. Our conclusions rely on symmetry 
reasoning, and only numerical details depend on the choice of the crystal field model which 
we use to illustrate the general arguments. 

The identification of the low-pressure, low-temperature hidden order of URu2Si2 is of 
basic interest. Starting from the high-temperature tetragonal phase, a symmetry-breaking 
transition can lead to an ordered phase with the following choices for the local order pa- 
rameter: A 2u and E u dipoles, Bi g , B 2g , and E g quadrupoles, B iu , and B 2u octupoles, a A 2g 
hexadecapole, and an Ai u triakontadipole 4j. 

It was always clear that the primary order parameter of URu2Si2 cannot be dipolar. The 
possibility of quadrupolar ordering has been extensively discussed Higher multipoles 
have been mentioned in a general context but have not been studied in detail. 

A recent /iSR study Q] finds that the symmetry of hidden order is different from A 2u (Q,) 
which is the symmetry of the high-pressure antiferromagnetic phase (the same structure 
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was ascribed to the supposed " micromagnetism" of URu2Si2, which is now understood to 
be extrinsic). The present experimental status is that the intrinsic low-pressure behavior of 
URu2Si2 is purely non-magnetic. Furthermore, a number of recent experiments proves that 
the hidden order breaks time reversal invariance, so it cannot be quadrupolar [7, 20, 3^|. In 



particular, Yokoyama et al |7[ found that uniaxial stress (which is time reversal invariant) 
induces large-amplitude antiferromagnetism, which breaks time reversal invariance. It is 
clear that stress must have acted on a medium which itself was non-invariant under time 
reversal: it must have been the octupolar phase j^]. 

We emphasize that stress-induced antiferromagnetism arises only if the stress is uniaxial, 
and perpendicular to the (001) direction. For tetragonal symmetry, octupoles (B iu and 
B 2u )i and dipoles (A 2u and E u ) are of different symmetry and therefore they do not mix. 
It follows that hydrostatic pressure cannot induce antiferromagnetism unless the pressure is 
high enough to .ead to a contp.etely different (pure.y dipo.ar) phase via a first-order phase 
transition. This was found in Ref. p. In contrast, uniaxial pressure perpendicular to the 
tetragonal main axis lowers the symmetry to orthorombic, allowing the mixing of dipoles 
and octupoles. 

We postulated that the hidden order is Tjr staggered octupolar order (Sec. |HJ). Uniaxial 
pressure cr|| (100) leads to the appearance of J z dipolar order of the same periodicity. The 
model works the same way if we postulate T xyz staggered octupole order, in which case a 
stress er||(110) gives J z antiferromagnetism. Since the octupoles Tjr and T xyz belong to dif- 
ferent one- dimensional irreps of the tetragonal symmetry {B 2u and B\ Ul respectively), in our 
theory a homogeneous system can show only one of the stress-induced effects. We hypothe- 
sized that the observed near-equivalence of the stress effect in (100) and (110) directions [7| 
reflects the presence of both kinds of order in a multi-domain structure. 

The same assumption about Tjr octupolar order explains the behavior in applied magnetic 
field (Sec. IHIj) . A field B||(001) mixes Tf octupoles with 0\ quadrupoles. Symmetry 
breaking is well-defined in the presence of magnetic field, and the transition to hidden order 
(now a mixed octupolar-quadrupolar order) remains second order up to a critical field _B cri i 
where T (B) — ► 0. 

We illustrated the symmetry arguments on the example of a crystal- field model (Sec. lIVp . 
The model has two versions: low-energy phenomena can be described by using three low- 
lying singlets, while for high energies (or fields, or temperatures) we need five states (the 
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previous three singlets plus a doublet). The three singlets are the same as in Santini's work, 
but their sequence was chosen to give an octupolar matrix element between the ground state 
and the first excited state. The presence of the doublet level is not essential at low fields 
(and low temperatures) but it splits in a magnetic field B||(001), and for a range of high 
fields, even weak quadrupolar coupling can give quadrupolar order which competes with the 
low-field octupolar order. We argued that the high-field order observed between 35T and 
38T 0Q 

is of quadrupolar nature, with a symmetry different from that of the low-field 

order. 

To conclude, we presented arguments showing that octupolar order of either B 2u or B\ u 
symmetry is the zero-field "hidden order" of URu2Si2 at ambient conditions. We limited 
the discussion to strictly on-site order parameters in a localized electron model with stable 
5/ 2 valence. However, within this restriction our scenario is compatible with the present 
knowledge about the phase diagram in the temperature-pressure-field space. The time 
reversal invariance breaking nature of the order is manifest in the effect that uniaxial pressure 
applied in certain directions can induce large-amplitude antiferromagnetism. 
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APPENDIX: ORTHOROMBIC SYMMETRY 



Quadrupolar moments couple to external stress. For instance, applying stress <r||(100), 
0\ quadrupolar moments are induced, while cr||(110) induces O xy . At the same time, the 
application of uniaxial stress lowers the symmetry from the tetragonal to one of its 
subgroups, changing the symmetry classification of all order parameters. This effect is 
described below. 

When we apply uniaxial pressure in direction (100), the previous C4, 5*4, C 2 ' and a v 
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cease to be symmetry operations and the residual symmetry is described by the group V 2 h- 
The corresponding classification of the order parameters is given in Table IIV1 We observe 
that under the new symmetry, the Tjr octupolar and the J z dipolar moments mix with 
each other, and this means that if the system possesses spontaneous staggered 7~r octupolar 
order, applying cr||(100) stress induces staggered J z dipolar moments. 

TABLE IV: Symmetry classification of the local order parameters for cr || (100) (left part) and 
for cr || (110) (right part). The subscripts g and u mean the even and odd under time reversal. 
H° = J x J y (J x — Jy) is one of the nine hexadecapoles. 
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When the uniaxial pressure is applied in (110) directions, it induces O xy quadrupoles. 
Now C4, 5*4, C 2 and ad have to be omitted from the symmetry group which is again the D 2 h 
point group, only comprised of different elements than in the cr || (100) case. For the present 
cr 1 1 (110) case, the symmetry classification of the order parameters is shown in Table HVl Now 
the T xyz octupolar moment mixes with the J z dipolar moment (if T xyz is staggered, so is J z ). 

When we apply uniaxial pressure along the tetragonal main axis (001), there is no sym- 
metry reduction, the original symmetry classification of the order parameters (Table H} 
remains valid. Neither the T xyz nor the Tf octupolar moments can induce J z magnetic 
moments, since they all correspond to different irreducible representations of the point 
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group. An analogous statement holds for the staggered moments. 
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